Spontaneously localized photonic modes due to disorder in the dielectric constant 
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We present the first experimental evidence for tlie existence of strongly localized photonic modes 
due to random two dimensional fluctuations in the dielectric constant. In one direction, the modes 
are trapped by ordered Bragg reflecting mirrors of a planar, one wavelength long, microcavity. In the 
cavity plane, they are localized by disorder, which is due to randomness in the position, composition 
and sizes of quantum dots located in the anti-node of the cavity. We extend the theory of disorder 
induced strong localization of electron states to optical modes and obtain quantitative agreement 
with the main experimental observations. 

PACS numbers: 78.67.Hc., 42.55.Dd, 78.55.Cr, 68.37.Uv. 



I. INTRODUCTION 

Localization of a photonic mode by disorder in a sta- 
tistically uniform and isotropic medium has been stud- 
ied for more than two decadesi"— , both in the optical 
and microwave spectral ranges. In this paper we pro- 
pose and demonstrate a new mechanism for disorder in- 
duced strong three dimensional (3D) localization of light. 
We achieved this in planar microcavities (PMCs), con- 
taining a layer of strain-induced self- assembled quantum 
dots (QDs)^"— . The dispersion of light in the PMC plane 
is modified by the distributed Bragg reflectors (DBRs) 
that trap the modes in perpendicular direction to the 
cavity plane. The QDs constitute an active material 
which, under proper excitation, emits photons, due to 
recombination of QD confined excitons^ii^. Randomness 
in the QDs in-plane position, composition and sizes pro- 
duces fluctuations in the system's dielectric constant. 
This randomness provides the necessary "attractive po- 
tential" , which, together with the modified dispersion, 
brings about strong photon localization in the PMC 
plane. Our novel mechanism differs from previous stud- 
ies, since it provides genuine disorder-induced localiza- 
tion rather than leaky (resonant) modesi"— . It also dif- 
fers from the "transverse localization"—, where the wave 
freely propagates in one direction while being trapped, 
by disorder, in the transverse directior^i^. We directly 
demonstrate the 3D localization of the photonic modes 
by measuring their photoluminescence (PL) intensity dis- 
tribution. 

It is worthwhile to emphasize from the start that the 
comparatively small size of the excitation spot (30 /im, 
or smaller) plays a crucial role in the interpretation of the 
experimental data. The point is that in the absence of 
disorder the shape of the PL spectrum emerging from the 
cavity is quite sensitive to the size of the excitation spot. 
If the size is sufficiently small, the PL spectrum becomes 



very narrow, in marked difference with the broad shoul- 
der predicted by the standard 'basic cavity physics'. If 
there had been no disorder, the experimental conditions 
of our paper would have resulted in an extremely narrow 
PL peak. The significant broadening of the spectrum, as 
observed in our experiment, is due entirely to disorder. 
The broadening occurs both towards frequencies below 
and above the cavity frequency (localized and resonant 
photonic modes, respectively). All this will be discussed 
in detail in the theoretical part of the paper. 



II. EXPERIMENTAL 

Our sample was grown by molecular beam epitaxy on a 
[1,0,0] oriented GaAs substrate. The sample consists of a 
PMC formed by a GaAs one- wavelength resonator, sand- 
wiched between DBRs made from 23(26) top (bottom) al- 
ternating GaAs and AlAs quarter- wavelength-layer pairs. 
As an active material in the resonator's anti-node, a sin- 
gle strain-induced self-assembled Ino.3Gao.7As QDs layer 
with average density of 5x10^ QDs per cm^ is usecU^. 
The use of such composition results in large and asym- 
metric QDs with typical lengths and widths of about 100 
and 30 nm respectively, as can be seen in the inset to 
Fig. [TJi. 

We used two methods to study the PMC. a) a diffrac- 
tion limited confocal optical scanning microscopy ( a x 60 
objective with numerical aperture of 0.85, obtaining spa- 
tial resolution of ~ 0.7 fim). b) a near field scanning 
optical microscopy (NSOM- Nanonics Cryoview 2000^^^ 
with spatial resolution of ^ 0.25//m). By applying either 
method we measure the lateral distribution of the elec- 
tromagnetic field above the sample surface. Two modes 
of excitation were used. In the first, the excitation was 
focused to a diffraction limit by the collecting objective 
— . In the second (the only one used for NSOM), the 
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excitation was focused at an oblique angle to a spot di- 
ameter of about 30 fim. Due to the high energy used for 
the excitation (632.8 nm light of a HcNe laser) the im- 
ages produced by the two excitation modes were almost 
identical. In Fig. [T^ we present PL spectrum (upper, 
black solid line) from a single point on the surface of the 
PMC sample together with a PL spectrum which was 
taken from the same sample (lower, solid blue line) after 
the upper DBR mirror was completely etched. The spec- 
trally broad blue PL line is centered at 1.33 eV and has 
full width at half maximum (FWHM) of approximately 
10 meV. This spectral line results from s-shell electron- 
hole recombination within the inhomogeneous population 
of QDs of various sizes and compositions. The radiative 
linewidth of a QD is ^ 2^eV— . 
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FIG. 1: a) PL spectrum of the PMC with (solid black 
line) and without (solid blue line) the upper DBR. The 
left inset represents spatially integrated PL spectrum from 
a 135^mx 200^771 area on the sample surface. b)-d) selective 
wavelength images of the area for 3 different spectral domains 
indicated by the green zones in a). The right inset is a scan- 
ning electron micrograph of the QDs layer. 

The rather symmetrical shape of the spectral line in- 
dicates that charge carriers are not distributed thermally 
between the QDs. Therefore, it is quite safe to assume 
that the spectrum accurately represents the actual en- 
ergy distribution of the emission from individual QDs. 
The spectrally sharp PL line in the figure is centered at 
1.319 eV and has FWHM of 80 jueV. As we explain below, 
in a disorder free PMC this line width should have been 
about an order of magnitude narrower due to the finite 
excitation area and diffraction limited collection spot. 

The asymmetric spectral shape of the line and its ac- 
tual linewidth which represents a Q-factor of ^15000, in 
similarity to earlier measurements^"—, are due to the dis- 
order in the dielectric constant as we quantitatively show 
below. 

In the left inset to Fig. [T^ we present on a semi loga- 
rithmic scale a spatially integrated ('far field') PL spec- 
trum from a 135^m wide x200/im long area on the sam- 
ple's surface. The spectrum was obtained by summing 
up large number of individual PL spectra, each obtained 



from a diffraction limited areal spot. The dominating 
spectral feature is a large, asymmetrically broadened line 
around the energy of the PMC mode. The line decays 
exponentially towards lower energies with a characteris- 
tic energy of 150 ± 30/iey. Towards higher energies the 
line decays with a characteristic energy of 380 ± 30fieV. 
The ratio between these characteristic energies agrees 
well with the prediction of our theoretical model for the 
disorder in the dielectric constant (see below). 

In Figs. [TJj-lTJi selective wavelength images of the 
scanned area are presented. The spectral domains which 
are used for each image are indicated by the green zones 
imposed on the 'far field' spectrum in Fig. [T^. The im- 
age in Fig. [TJj is obtained within a 0.5 meV energy win- 
dow containing the PMC mode. It shows almost evenly 
distributed emission from the surface. In contrast, the 
image in Fig. [1];, which is obtained within a 0.5 meV 
window, 1 meV below the PMC mode, shows emission 
emerging from randomly distributed spots on the sur- 
face. Likewise, the image in Fig. [TJl, which is obtained 
from energy window of same width, located 3 meV lower, 
shows only 2 bright centers, from which all the emission 
results. 
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FIG. 2: a)-e) middle panels - selective wavelength NSOM 
images of an area in the vicinity of the lowest energy mode 
observed in Fig. [T] The spectra (left panels) are obtained 
from the brightest area pixels indicated by the + signs on the 
images. The green areas on the spectra, indicate the spectral 
domains used for generating the images. The right panels 
represent analytically calculated ID intensity distributions of 
the confined optical modes obtained by solving Eq. (I2J with 
the parabolic dielectric constant fluctuation model presented 
in the right panel of e) where Ae/e = 0.7% and e the dielectric 
constant of GaAs. The calculated mode energies are presented 
by dashed vertical lines on the spectra to the left. 



In Fig. [3] we turn our attention to one of the low energy 
electromagnetic modes seen in Fig. [TJi. Such localizing 
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centers were observed in various mode energies and var- 
ious locations on the sample surface. In the right panels 
of Figs.[2K-l2l3, we display selective wavelength images ob- 
tained by scanning the Aluminum coated fiber tip of our 
NSOM microscope in tapping mode over the sample sur- 
face. The corresponding spectra in the left panels, are 
obtained from the brightest area pixels indicated by the 
+ signs on the images. The green areas on the spectra 
indicate the spectral domains used for generating these 
images. These spectra and corresponding images reveal 
the following observations: Four different localized pho- 
tonic modes are observed at energies below that of the 
PMC mode. Each mode is composed of two cross lin- 
early polarized equally intense components (not shown). 
The components are ^ 80^6^ apart and the lower en- 
ergy one is polarized along the [0,-1,1] dircctioniS. The 
lowest energy mode is about 5 meV below the energy of 
the PMC mode. The in-plane spatial intensity distribu- 
tion of this mode, as depicted by the image in Fig. 
is node-less, and has an elongated shape oriented along 
the [0,-1,1] crystallographic direction, along which the 
QDs are also elongated (see inset to Fig. [1]). The mode 
dimensions are roughly 3x2fim^. The peak intensity is 
about 40 times stronger than that of the PMC mode. The 
next three images in Fig. [5Jd, c and d, present elongated 
spatial intensity distributions composed of one, two and 
three nodes, respectively. Each mode's energy is ~1.5-2 
meV higher than the energy of the preceding mode. The 
selective wavelength image in Fig. ^ is obtained at an 
energy window, which resides on the higher energy side of 
the PMC mode. It displays evenly distributed emission 
from everywhere except from the actual locus of the lo- 
calized modes, where extended modes are forbidden. The 
magnitude and shape of the spot from which emission at 
this wavelength is missing is an estimate for the spatial 
extent of the region which localizes the optical modes. 
The region is comparable in shape and size to the most 
localized, lowest energy mode. The localized modes are 
insensitive to temperature and excitation density varia- 
tions, hence they do not result from thermal or non-linear 
effects. The emission is linear with the excitation den- 
sity between 1-10^ Watt/cm^. Increasing the tempera- 
ture from 10 to 80K results in a small, rigid, red spectral 
shift, as expected from the temperature dependence of 
the material bandgap. 



III. THEORY 

Since the sizes of the spots in Fig. 2b-2d significantly 
exceed the size of an individual QD, the only way to 
account for the experimental observations is to assume 
that each spot corresponds to a localized electromagnetic 
mode. Therefore the low frequency tail of the PL spec- 
trum in Fig. la should be interpreted as a tail of localized 
modes with energies below that of the PMC mode, Wc, 
similar to the tail of localized electron states formed be- 
low the band-edge due to disorde r^^i^^ . To demonstrate 



how the presence of a cavity leads to the in-plane trap- 
ping, we consider the scalar wave equation corresponding 
to the frequency fl: 



£o + Ss{p,z) 



* = 0, 



(1) 



where £o is a uniform dielectric constant inside the cav- 
ity, p = {x,y), and Ss(p,z) is the perturbation due to 
the randomly positioned QDs that reside in a narrow 
layer, of width (5z ~ 5 nm, in the middle of the cav- 
ity. Let us emphasize that an individual QD is much too 
small to be able to trap a photon. Actually, in the con- 
tinuum description employed in Eq. ([T|) individual QDs 
are not resolved. What matters are the fiuctuations in 
the areal density of the QDs, which manifest themselves 
in fiuctuations of the dielectric constant Se. Since the 
electromagnetic modes are confined in the z-direction to 
a much shorter length than in the plane, the function 
\I'(p, z) can be factorized as $(p)cos(7rz/(i), where d is 
the cavity width, and $(p) satisfies the equation 



(2) 



In this equation, 6e{p) is the effective in-plane fluctua- 
tion of the dielectric constant obtained from Se{p, z) by 
averaging over the z-distribution of the cavity mode in- 
tensity, namely, 5e{p) = | /^^yj ^^iP^ cos^(7rz/d). The 
analogy between Eq. ^ and the Schrodinger equation 
makes it clear that, for positive 5e{p) , Eq. ([2]) admits 
localized solutions, with frequency VL smaller than the 
cavity frequency Wc = l-Kcj d^fe^. Thus, even a small 
enhancement in Se{p) acts as an attractive potential lo- 
calizing the mode laterally. To substantiate that this 
mechanism of trapping is relevant to our experimental 
observations, we note that the spatial intensity distribu- 
tions and the energies of the modes in Fig. [5] are quite 
accurately fitted by solving Eq. ([2]) using a one dimen- 
sional parabolic model (^Ss{x) = Ae[l — {ax)'^], \ax\ < ij 
describing a local fluctuation in the dielectric constant 
(see Fig. [2j;, right panel, where the best fltted values are 
a — 0.3^m~^ and Ae/eo = 0.7%). This agreement, and 
the fact that even in the case of the most localized modes, 
differently polarized components have almost the same 
energy, justifles the use of the scalar approximation. 

To flnd the shape of the tail of the localized modes 
distribution, we have to specify the disorder "potential" , 
5e{p). The simplest possible model is the "white-noise 
potential" which amounts to treating the QDs as point- 
like resonant units, randomly distributed in the plane 
(this simple model does not describe the anisotropy ob- 
served in the inset of Fig. [T] as discussed below). Fluctu- 
ations in the local dielectric constant are given as 



5e{p) ^aJ 



(3) 



where SNi^{p) is the deviation of the number of QDs 
(per unit area near p), with resonance frequency in the 
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interval du, from the average value NPu{u>)duj. Here N 
is the total areal density of the QDs, regardless of their 
frequency and Pd{'-^) is the probability distribution of 
the resonance frequencies (the solid, blue line in Fig. [TJ. 

To estimate the constant A, we consider the f2 — > 
limit and imagine that the entire layer consists of QDs, 
i.e. of Ino.3Gao.7As. In this case, the enhancement of 
the dielectric constant in the layer, as compared to the 
surrounding material of GaAs, is (e/n^s — sguAs) ■ 0.3. 
This value should be multiplied hy 5 = 5/300, which is 
the ratio of the QDs layer width (5 nm) to the cavity 
width (^300 nm). There is an extra factor of 2, since 
the QDs are located at the anti-node of the PMC mode. 
Thus, the effective enhancement, due to the presence of 
InAs, is Ae = 2.25 x 10~^. The constant A can now be 
estimated from Eq. ^ a.s A = 2.25 x 10~^S'a;g, where S 
is the typical area of a QD and wo ~ Wc is the central 
frequency of the distribution Pd{u!). Note that the cm- 
ployed averaging in the z-dircction, as well as our use of 
the scalar approximation is justified by the smallness of 
the parameter S. Indeed, the effect of disorder-induced 
coupling between different cavity modes is proportional 
to 5^, whereas the effect of the QD layer on a given mode 
is proportional to S. 

The function 5N^{p) is a complicated, random func- 
tion of position p and frequency uj. Using the white- noise 
model for 5N^ (p) and Eq. ([3]) , we obtain the probability 
P {6e{p)} for a particular realization 6e{p) of the fluctu- 
ating part of the dielectric constant 



P {ds (p)} '-^ exp ■ 



2l^2 



P£,(yj)dio 
NA^ \J {uj - f})2 



(4) 

Formally, the integral over lo in Eq. ^ diverges at 
UJ il. In this regard, it should be noted that the inte- 
gral over u! is actually an approximation for a sum over 
discrete values, Wi, of the resonant frequency of the dots 
in the relevent area (of order 1 /im^). The meaning of 
this area can be traced back to Eq. ([3]) which defines a 
macroscopic (albeit still fluctuating in space) dielectric 
constant. The spatial scale, over which 5e{p) is chang- 
ing, should be at least few times larger than the "micro- 
scopic" distances (size of the dots and their separation, 
of the order of 0.1 fxm). Thus, the minimal scale at which 
Se{p) can be meaningfully defined should be in the sub 
micron range. We are interested in f2 below Wc, i.e. in 
frequency which belongs to the tail of the distribution 
Pd(oj)- a simple estimate shows that for such VL the 
probability to encounter a dot with LOi close to fi is ex- 
tremely small. Moreover, for such f2 the integral over lo 
in Eq. (jH) can be, roughly, replaced by (ojo — VL)^'^ (note 
that such a replacement would not be possible for f2 in 
the "bulk" of the probability distribution Po{(jj))- The 
two equations, ^ and define the statistical prob- 
lem of finding the probability that a localized mode, in 
a given frequency interval, will be created. This problem 
is completely identical to the two-dimensional versionii 
of the problem of the tails of electron states in a ran- 



dom potentiali^ii^. The final result for the probability 
V{il)dQ,dS to find a trap, in an area dS, which can sup- 
port a localized mode in the frequency range dfi can be 
presented as 



£o^c 
5- exp 



, (5) 



where the pre-exponcntial factor corresponds to the den- 
sity of the PMC modes at frequencies > ujc- Eq. ([5]) 
can be viewed as a qualitative interpolation between the 
maximal value of the modes density at = Wc, and the 
low-density "tail". A more detailed crossover behavior 
between the tail at i7 < Wc and V = EquJc/'^c^ = Pa at 
n > (jjc can be found in Rcf. [l^ Important is that this 
crossover is determined by a single scale of frequencies 



A^N 



UJt = i^c 



327reoc2(wo — (^c 



(6) 



which follows from Eq. ([5|), so that in the tail region 
V{fl) can be rewritten as V{^1) « Va exp [— (wc — ^)/^t]- 
In Eq. ([6|) we replaced {ojq — il) by {ujq — Wc), since, 
experimentally, the difference (wq — Wc) is much bigger 
than the tail width, wt. 

With the above estimate for the constant A, Eq. ([5]) 
yields a characteristic energy of lOOztTOfieV for the initial 
decay of the PL intensity below HuJc- It agrees well with 
the measured value (dashed red line in the left inset to 
Fig. [T^). The relatively large uncertainty results from 

I the uncertainties in the QDs density and their average 

^area. 

' Up to now we have focused on explaining the PL spec- 
trum tail for id < ujc- In this domain, the PL intensity, 
/(rj), outside the cavity is simply proportional to the 
density of trapped modes, V{^)- In other words, the Q- 
dependence of the light intensities inside and outside the 
cavity is the same for uj < uJc- This is not the case for the 
domain of frequencies uj > ujc to which we now turn. To 
calculate in this domain, we start by recapitulating 
the basic microcavity physics, in the absence of disorder. 
In this case the emission spectrum is controlled by the 
Airy factor— 



Ai{n,e) = 



T 



1 — fir2 exp [2i(/)(ri, ( 



(7) 



where 9 is the angle at which the radiation exits the cav- 
ity, f2 > is the radiation frequency, 7'i (r2) is the re- 
flection amplitude from the upper (lower) Bragg mirror 
and r = 1 — is the transmission coefficient (ri and r2 
are taken to be real). The phase (f>{i^, 0) = ^Jeo^dcosO/c 
is acquired in course of propagation between the bound- 
aries z = ■izd/2. At frequencies, f2, close to uJc Eq. ([7]) 
has a sharp maximum aX 6 ^ 6^, where the angle On is 
defined as 



2{n - UJc 



1/2 



(8) 
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so that Ai{n,9) can be simplified to 



Ai{n,e) = 



T 



T2 + 47r2( 



(9) 



Note that condition = 9q coincides with the dispersion 
law of the waveguided mode. This follows from Eq. ^ in 
the absence of disorder. Indeed, for in-plane wave vector, 
q = y/eo^siii9/c, Eq. ([2]) yields q w y/2eouJc{fl — lu^^/c, 
leading to 9 = 9a for 9^ ^ 1. 

It follows from Eq. (j9l) that the emission spectrum, 
/(f2), has a form of a plateau, which corresponds to the 
basic microcavity physics. Indeed, integration over 9 
yields 



CJO 

oc J d9 9Ai{n,9) = i 



2 /27r6l2 

- arctan 

V T 



TT 



(10) 

The second term describes a slight smearing of the 
plateau near fl = ujc- 

Let us consider now photons which exit at an angle 
9, from some small collection spot. It is important to 
realize that there is a correlation between the angle 9 
and the average distance L{9) (measured from the col- 
lection spot) at which the corresponding photons have 
been created. Indeed, a typical photon wih bounce 1/T 
times before escaping the cavity. Therefore photons that 
bounce at angle 9 and escape at the collection spot must 
be created at a distance 



9rl 

L{9)^T-^dian9K—. 



(11) 



Thus, the finite size of an excitation area can serve as an 
efficient cutoff for a maximal escape angle 9m ■ Assuming 
a circular excitation area, of radius R, and choosing a 
small collection spot at the center of the circle, one finds 
a maximal escape angle 



RT 



(12) 



Of course, this relation is valid only if (and there- 
fore R) is not too large, namely, it should be smaller 
than the maximal escape angle for a fully excited cavity. 
Since the angular cutoff translates into a frequency cutoff 
^Ira = uJc[^ + \9l^) (see Eq. ([8])), we arrive at the impor- 
tant conclusion that the frequency range of the collected 
photons can become very narrow, if the excitation area 
is small enough. In our experiment T = 0.44 x 10~^, 
R = 30/x, so that the collected radiation constitutes a 
narrow peak at near lJc, of relative width 




1 X 10" 



(13) 



Obviously, this peak can be replaced by a (5-function 

iin)^ioSin~iOc). (14) 



Qualitatively, the meaning of Eq. ([14)) is that, due to the 
finite excitation spot, only the photons emitted normally 
to the mirror emerge outside the cavity. Quantitatively 
we note that the relative width of the experimentally 
measured PMC mode (« 1 x 10^"*) is about an order of 
magnitude wider than the one calculated by Eq. (|T5)) . As 
we show below, this broadening and the actual spectral 
shape of the PMC mode, results from the disorder. The 
replacement of the step-like spectrum of the clean cavity 
(for infinite excitation spot) by the i5-function Eq. 
(for finite excitation spot) does not affect the emission 
spectrum of disordered cavity in the domain fl < ujc- In 
this domain the dependence, /(ri), is still determined by 
Eq. (O. However, in the domain > Wc, the (5-peak 
for the clean cavity gets broadened due to scattering of 
extended modes by the disorder. 

In the presence of disorder, the above reasoning, which 
is based on the ray picture, is not applicable since all 
the solutions, $^(p), of Eq. ^ corresponding to the 
eigenstates, Vt^, have components with q = (in the ray 
picture they correspond to rays with normal incidence to 
the microcavity plane). The amplitude of this component 
is is given by jd^p $p(p). Therefore, the shape of the 
spectrum acquires the form /(fi) = IqA{^) where the 
function A{Vl) is defined as 



AiSl) 



J2 ld^p\<i>M^s{n-n,). 



(15) 



We note that the function A{^1) as defined by Eq. ([T 
is normalized 



dnA{n) = 1, 



(16) 



Thus, the disorder broadens the S-peak (Eq. ([T^ ) with- 
out changing its area. Note also, that the calculation 
of the spectral shape, A{fl), is formally reduced to the 
calculation of the spectrum of two-dimensional exciton 
absorption in the presence of a white-noise random po- 
tential. This calculation was done earlier in Ref. |20| . 
where the zero-momentum components of the exciton 
states were selected by the matrix element for light ab- 
sorption^S. 

It is also important to note that the spectral shape, 
A{il), is a function of the ratio (fi — uJc)/uJt, where Wt is 
the same characteristic frequency, (Eq. ([5])), which de- 
termines the spectral shape of the tail of the localized 
photon states. Indeed, deep in the tail, (ujc — fi) 3> wt, 
all the trapped solutions, ^^{p), at a given Vl^ — ^ 
are identical to each other, and Eq. (|15p is essentially 
reduced to the density equation of the trapped modes 
(xexp[-(a;c-0)/wt] (Eq. (O). 

In the opposite limit (O — Uc) 3> the density of 
photon states is constant. Here, therefore, the decay 
of A{^) is determined by the fall-off of the integral 
jd^p (f>^(p). This integral can be estimated to first or- 
der by perturbation expansion in 6e{p). Thus it is pro- 
portional to {n — uJc)~^ ■ Accurate calculation^ yields 
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^(fi) w — uJcY . This way, the width of the emis- 

sion spectrum was found to be S.Twj. 

The overall spectral shape of A{^) was found in Rcf. 
[20by combining the normalized (Eq. ([IGp ) low- frequency 
and high-frequency asymptotes. The calculated full 
width at half-maximum (FWHM) of the emission spec- 
trum is S.Twt, in excellent agreement with the experi- 
ment, where it is measured to be S.Swt. The calculated 
ratio between the high energy and low energy slopes is 
2:1 in very good agreement with the measured ratio of 
2.5. 



IV. DISCUSSION 

While the white noise model is successful in explaining 
the initial drop of the density of localized modes, it can- 
not account for the experimental data in the deeper tail, 
starting around O = 1.318 eV and below. It does not ex- 
plain the non-monotonic decay of the probability as seen 
in Fig. [1^. It fails in predicting the observed multiplicity 
of confined modes which belong to the same localization 
region (Fig. [2]) , and it does not explain the anisotropy in 
the electromagnetic mode shapes. 

The failure of the white noise model, for O < Wc, is 
not surprising. Indeed, when the mode frequency be- 
comes smaller, localization of the mode becomes tighter. 
Thus, various local features of the random potential - 
such as correlations in Se{p), the "granular" structure of 
the disorder and the inherent anisotropy of the system - 
become more important. Obviously, all these features are 
not captured by the entirely featureless, universal white 
noise model. 

Proper account of these specific features of randomness 
might explain some of the data for the deeper H-tail. 
First, once correlations are introduced, the most prob- 
able deep fluctuation is no longer a fluctuation which 
captures one mode only. Second, the probability of find- 
ing a localized mode common to two localizing centers, 
rather than to one, may prevail at a given energy below 
the energy of the PMC mode. Indeed, localized modes 
between 1.317 eV and 1.318 eV always exhibit a spa- 
tial structure consisting of two bright spots separated 
by a dark region. This observation suggests existence of 
traps of a double- well shape. The origin of such traps, as 
well as the multimode traps in the deep tail (below 1.316 
eV), may be traced to the granular nature of the disor- 



der. Although, on the average, the QDs are uniformly 
distributed, there exist configurations with several dots 
coming close to each other, within a distance of a dot size. 
Such rare configurations constitute more efficient traps 
than the white-noise fluctuations and therefore dominate 
in the deeper tail. The clusters of QDs arc akin to clus- 
ters of potential wells in the Lifshitz model of a disor- 
dered electronic system3ii^. Thus, anisotropic QD clus- 
tering, can account for the shape and multiplicity of the 
strongly localized modes in the deeper tail. We are un- 
able to make quantitative predictions for the density and 
spatial structure of these modes, since no comprehensive 
theory exists. Yet, we show that dielectric constant en- 
hancement of order 1% over distance of order Ifim, does 
account for the experimental observations (Fig. . 

In conclusion, we demonstrate disorder-induced 
trapped photonic modes in a microcavity with an em- 
bedded layer of QDs. In the lateral in-plane direction the 
modes are localized by spatial fluctuations of the dielec- 
tric constant due to randomness in the location and com- 
position of the QDs. Our theory emphasizes the universal 
features of the trapping mechanism: the necessary combi- 
nation of mirrors in the vertical direction and the lateral 
disorder. The theory quantitatively estimates the char- 
acteristic energy by which the PL intensity decays below 
the PMC mode. Since the non universal, system-spcciflc 
details of the randomness are not accurately known, the 
spatial shapes and energies of the modes deeper in the 
energy tail, are qualitatively discussed only. 

In addition we show both experimentally and theo- 
retically that if the size of the excitation spot is suffi- 
ciently small, as it is in confocal and near fleld scanning 
microscopy, the microcavity emission spectrum becomes 
very narrow, in marked difference with the prediction of 
the standard cavity considerations, which apply to uni- 
form excitation. The disorder broadens this spectrum, as 
well. Thus, we quantitatively measure and calculate the 
disorder induced broadening below and above the cavity 
mode frequency. 
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